Structure of Fermionic Density Matrices: Complete iV-representability Conditions 
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We present a constructive solution to the iV-representability problem — a full characterization of 
the conditions for constraining the two-electron reduced density matrix (2-RDM) to represent an 
iV-electron density matrix. Previously known conditions, while rigorous, were incomplete. Here we 
derive a hierarchy of constraints built upon (i) the bipolar theorem and (ii) tensor decompositions 
of model Hamiltonians. Existing conditions D, Q, G, Tl, and T2, known classical conditions, 
and new conditions appear naturally. Subsets of the conditions are amenable to polynomial-time 
computations of strongly correlated systems. 
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PACS numbers: 31.10.+Z 

The wavefunction of a many-electron quantum sys- 
tem contains significantly more information than neces- 
sary for the calculation of energies and properties. In 
1955 Mayer proposed in Physical Review computing the 
ground-state energy variationally as a functional of the 
two-electron reduced density matrix (2-RDM) which, un- 
like the wavefunction, scales polynomially with the num- 
ber N of electrons [l|-|3j]. However, the 2-electron den- 
sity matrix must be constrained to represent a many- 
electron (or _/V-electron) density matrix (or wavefunc- 
tion); otherwise, the minimized energy is unphysically 
below the ground-state energy for N > 2. Coleman 
called these constraints N -representability conditions 
and the search for them became known as the N- 
representability problem [Hj-Tiol) . In 1995 the National 
Research Council ranked the iV-representability problem 
as one of the top unsolved theoretical problems in chem- 
ical physics [11]. While progress was limited for m any 
years, recent advances in theory and optimization 12H20l| 
have enabled the application of the variational 2-RDM 
method to studying strong correlation in quantum phase 
transitio ns 12111 . quantum dots 22j, polyaromatic hydro- 
carbons [23|, firefly bioluminescence 24 1, and metal-to- 
insulator transitions [25j | . 

Despite the recent computational results with 2-RDM 
methods, a complete set of 7V-representability condi- 
tions on the 2-RDM — not dependent upon higher-order 
RDMs — has remained unknown. While formal solutions 
of the N- representability problem were developed in the 
1960s [5j, |26j, practically they required the A^-electron 
density matrix [H, 0|. In this Letter we present a con- 
structive solution of the 7V-representability problem that 
generates a complete set of A^-representability conditions 
on the 2-RDM. The approach is applicable to generat- 
ing the iV-representability conditions on the p-RDM for 
& n y P < N. The conditions arise naturally as a hier- 
archy of constraints on the 2-RDM, which we label the 
(2, g)-positivity conditions, where the (2, 2)- and (2, 3)- 
positivity conditions include the already known D, Q, G, 
Tl, and T2 conditions 0, H, 0, E3 • The second number 
in (2,q) corresponds to the higher g-RDM which serves 



as the starting point for the derivation of the condition. 

A key advance in extending the (2, <7)-positivity con- 
ditions for q > 3 is the use of tensor decompositions in 
the model Hamiltonians that expose the boundary of the 
Af-representable 2-RDM set. The decompositions allow 
the terms in the model Hamiltonians to have no more 
than two-body interactions through the cancelation of all 
higher 3-to-q-body terms. A second important element 
is the recognition that when q = r where r is the rank 
of the one-electron basis set the positivity conditions are 
complete. The hierarchy of conditions can be thought of 
as a collection of model Hamiltonians @. For example, 
the 'basic' (2,2)-positivity conditions are both necessary 
and sufficient constraints for computing the ground-state 
energies of pairing model Hamiltonians [2, E3] , often em- 
ployed in describing long-range order and superconduc- 
tivity. 

Consider a quantum system composed of N fermions. 
A matrix is a fermionic density matrix if and only if it 
is: (i) Hcrmitian, (ii) normalized (fixed trace), (Hi) an- 
tisymmetric in the exchange of particles, and (iv) posi- 
tive semidefinite. A matrix is positive semidefinite if and 
only if its eigenvalues are nonnegative. The p-particle 
reduced density matrix (p-RDM) can be obtained from 
the iV-particle density matrix by integrating over all but 
the first p particles 
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The set of N D is a convex set which we denote as P N 
while the set P D is a convex set which we denote as Pfr, 
the set of A^-representable p-particle density matrices. A 
set is convex if and only if the convex combination of any 
two members of the set is also contained in the set 
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w) N D 2 €P N , 



(2) 



where < w < 1. The integration in Eq. ([1} defines 
a linear mapping from P N to P^, which preserves its 
convexity. 

The energy of a quantum system in a stationary state 
can be computed from the Hamiltonian traced against 
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the state's density matrix. For a system of N fermions 
we have 



E = Tr(H N D 



(3) 



If the Hamiltonian is a p-body operator, meaning that it 
has at most p-particle interactions, then the energy can 
be written as a functional of only the p-RDM 



E = Tr( H P D) 



(4) 



For a system of N electrons the Hamiltonian generally 
has at most pairwise interactions, and hence, the energy 
can be expressed as a linear functional of the 2-RDM. 
Except when N = 2, however, minimizing the energy as 
a functional of a two-electron density matrix 2 D G P 2 
yields an energy that is much too low. To obtain the 
correct ground-state energy, we must constrain the two- 
electron density matrix to be A-representable, that is 

2 D e P%. 

Based on the equivalence of the energy expectation 
values in Eqs. (J3]) and (|1]), we can use the set Pjy of N- 
representable p-particle density matrices to define a set 
Pjy of p-particle (Hamiltonian) operators p O that are 
positive scmidcfinitc in their trace with any A^-particle 
density matrix 



pP* = { p O\Tr{ p O p D) > for all P D G P^}. 



(5) 



The set Pjy is said to be the polar (or dual) of the set 
Pjy. Importantly, by the bipolar theorem 2(| 27 1, the set 
Pjy also fully defines its polar set Pjy as follows 

P p = { p D\Tt( p O p D) > for all p 6 G P^*}. (6) 



By Eq. ((6|) we have a complete characterization of the 
iV-representable p-R DMs from a knowledge of all oper- 
ators p O G Pjy* [26|. This analysis shows formally that 
there exists a solution to the ./V-representability prob- 
lem [j| [26|, but it does not provide a mechanism for 
characterizing the set Pjy*. 

To characterize Pjy*, we assume that the A^-fermion 
quantum system has r orbitals and hence, r — N holes. 
A convex set can be defined by the enumeration of its 
extreme elements, that is the elements (or members) that 
cannot be expressed by a convex combination of other 
elements 0, The definition of P p * in Eq. © for 

p < N can be extended in second quantization to include 
p > N 



P p * = { p O\Tr( p O D) > for all D} 



(7) 



with the p O being polynomials in creation and annihila- 
tion operators of degree 2p. Because in second quanti- 
zation the value of N is defined in the density matrices 
N D rather than in the operators p O (28|, the set Pjy* 
provides complete A^-representability conditions on the 




FIG. 1. The convex set P% of 2-body operators that are 
positive semidefinite in their trace with any iV-particle density 
matrix is contained within the convex set P%* of analogous 
3-body operators, which in turn is contained within the set 
Pn*. Hence, the extreme points of P% can be characterized 
completely by the convex combination of the extreme points 
of Pn* , which are given by Eq. (JSJ). 



p-RDM for any N between 2 and r. The extreme opera- 
tors in the set Pjy* can be written as Hermitian squares 
of operators (2i| 



r O,. = r d r C\ 



(8) 



where the r Ci are polynomials in the creation and anni- 
hilation operators of degree less than or equal to r (i.e., 
Eqs. and (f2T)]) b Because any operator P C with p > r 
reduces to a polynomial of degree r in its operation on 
any N D, the sets Pfj* with p > r do not contain addi- 
tional information about the positivity of N D. To estab- 
lish this reduction, we rearrange terms in P C of degree 
greater than r into a normal order with either more than 
annihilation operators to the right of the creation oper- 
ators or more than r—N creation operators to the right of 
the annihilation operators; in either situation, the terms 
of degree greater than r vanish in their operation upon 
any N D. 

The operators p O that constrain the p-RDM to be A^- 
representable in Eq. ([6]) are also necessary to constrain 
the g-RDM to be A^-representable where q > p; formally, 
each p O G Pjy* can be lifted by inserting the number 
operator to the (q — p) power to form a q O G P^* [14 1. 
Therefore, as illustrated in Fig. 1, we have the following 
set relations 



P 



N 



C P 3 * 



C P p * 



C P r * 



(9) 



Consequently, extreme operators r Oi of Pjy* can be com- 
bined convexly to produce all p-body operators p O G 
Pjy* , and hence, the extreme points of Pjy* can be char- 
acterized completely by the convex combination of the 
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extreme points of PL*. More generally, convex combina- 
tions of extreme q Oi £ P^* generate all p-body operators 
p O 6 P%f* for p < q. Depending upon the order of the 
creation and annihilation operators in r Oi, the normal- 
ordered terms will have either positive or negative coeffi- 
cients. Convex combinations of the r Oi can be chosen to 
cancel the coefficients of all terms of degree greater than 
p. Extreme elements are generated from the minimum 
number of convex combinations to effect the cancelation. 
This characterization of the set Pjy provides a construc- 
tive solution of the iV-representability problem for the 
p-KDM. 

The constructive solution — convex combinations of 
the operators in Eq. ([5]) — generates the existing N- 
representability conditions as well as new conditions. 
The (l,l)-positivity conditions Q are derivable from the 
subset of r Ci operators in Eq. ([§]) of degree 1 
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(10) 
(11) 



Keeping the trace of the corresponding one-body opera- 
tors 1 Od and 1 Oq against the 1-RDM nonnegative for 
all values of bj yields the conditions, 1 D >z and 1 Q >; 0, 
where 1 D and 1 Q are matrix representations of the 1- 
particle and the 1-hole RDMs and the symbol >; indicates 
that the matrix is constrained to be positive semidcfinite. 

Similarly, the (2,2)-positivity conditions [5| follow from 
considering the r Ci operators of degree 2 in Eq. flSJ) 



C D 


jk 


(12) 


Cq 


= 22 bjkO-j&k 
jk 


(13) 


C g 


= ^ bjkOjO, k . 


(14) 




jk 





Restricting the trace of the corresponding two-body op- 
erators 2 0_d, 2 Oq, and 2 Og against the 2- RDM to be 
nonnegative for all values of bj^ defines the conditions, 
2 D hO, 2 Q h 0, and 2 G y 0, which constrain the proba- 
bilities for finding two particles, two holes, and a particle- 
hole pair to be nonnegative, respectively. 



In general, the (q, g)-positivity conditions [12J, ll4| fol- 
low from restricting all g-body operators q O in Eq. © 
to be nonnegative in their trace against the q-RDM [14] . 
While the (q, q)-positive operators are not two-body op- 
erators for q > 2, convex combinations of them gen- 
erate two-body operators 2 <J 6 P% that enforce the 
iV-representability of the 2-RDM. We refer to necessary 
7V-representability conditions arising from convex combi- 
nations of (q, g)-positivity conditions as (2, g)-positivity 
conditions. 



The simplest such constraints, the (2,3)-positivity con- 
ditions, arise from keeping convex combinations of 3- 
body operators in Eq. ((SJ nonnegative; for example, 
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Ot\ — — (Cti.i Cti,i + ^TiaC T12 ) 



Ot2 — -(CT2,1 C T21 + Ct2,2C T22 ) 



where 



C 



Tl.l 



J2b jkl a]a{al 



Cti,2 = yb* M ajakai 
jki 

Ct2,i = bjkiajalai - 
jki 

Ct2,2 = Y b* ikl aja k a[ 
jki 



E^3 

3 



(15) 
(16) 

(17) 
(18) 
(19) 
(20) 



These conditions, known as the Tl and generalized T2 
conditions were developed by Erdahl Ml and implemented 
by Zhao at al. 15 1 and Mazziotti [14J. In general, they 
significantly improve the accuracy of the 2-positivity con- 
ditions. 

Although the constructive proof given above indicates 
that a complete set of iV-representability conditions can 
be generated from convex combinations of extreme ele- 
ments of fjy*, additional conditions have not been dis- 
covered beyond the (2,2)- and (2,3)-positivity conditions. 
For example, what about (2,4)-positivity conditions — 
that is, TV-representability constraints on the 2-RDM 
arising from convex combinations of 4-body operators in 
Eq. (|8|)? First, we derive a class of (3,4)-positivity con- 
ditions on the 3-RDM. 

Consider the nonnegativity of the following operator 
O formed by the convex combination of a pair of 4-body 
operators from Eq. §5§ 



O 



(Cxxxx C^,^ 



+ C'xoooC'xooo) 



(21) 



where the symbols x and o represent creation and annihi- 
lation operators, respectively, in the C operators defined 
as follows 



jklin 



(22) 
(23) 



jklra 



Importantly, the expectation value of O with dj 



klm 



Djklm 



requires the 4- RDM because the cumulant part A 



of the 4- RDM [l|, Ho| does not vanish 

,* (4,A3klm _4 A jqst s .. 
klmUpqst \ ^pqst ^pklm) r u - 



jklrnpqst 



(24) 
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TABLE I. A class of (2,4)-positivity conditions can be derived from convex combinations of the (4,4)-positivity conditions that 
cancel the 3- and 4-particle operators. We achieve the cancelation through tensor decomposition in the model Hamiltonians. 

(2,4)-Positivity Conditions 

Tr((3C'xxxxC'xxxx 4~ CxxxoCx X xo 4" C'xxoxC'xxox 4~ CxoxxC XO xx 4" CoxxxCqxxx 4" C'ooooC'oooo) ^ 

Tr((3C X xxoCx X xo 4~ CxxxxCxxxx 4~ CxxooC'xxoo 4" CxoxoCx OX o 4~ CoxxoCo X xo 4" C'oooxC'ooox) ^) ^ 

Tr((3C'xxoxC'xxox 4" C'xxoo^xxoo 4" CxxxxC X xxx + CxooxCxoox 4~ C OX oxCo X ox 4~ C'ooxo^ooxo) -^) ^ 

Tr((3C'xoxxC'xoxx 4* CxoxoC xoxo 4" C X ooxC xoox 4~ C X xxxCx X xx 4" C 00 xxCo OXX 4" C*oxoo^oxoo) ^) ^ 

Xr((3C'oxxxC'oxxx 4" C'oxxoC'oxxo 4" C'oxoxC'oxox 4~ C'ooxxC'ooxx 4~ CxxxxCxxxx 4~ C'xooo C'xooo ) D) ^ 

Tr((3C XX ooC xxoo 4~ C'xxoxC'xxox 4" C'xxxoC'xxxo 4" C XO ooC xooo 4~ C'oxooC'oxoo 4" C'ooxxC'ooxx) ^) ^ 

Tr((3C'xooxCxoox 4" CxoooCxooo 4~ C XO xxC xoxx 4" C'xxoxC'xxox 4~ CoooxCooox 4" CoxxoCoxxo) D) ^ 

Tr((3CxoxoCx OXO 4* CxoxxCxoxx 4" CxoooCxooo 4" CxxxoCxxxo 4" CooxoCooxo 4" CoxoxCoxox) ^) ^ 



To obtain additional iV-representability conditions re- 
quires that the dependence of the C operators on the ex- 
pansion coefficients be generalized from linear to nonlin- 
ear. Specifically, to obtain 3- RDM conditions beyond the 
(3,3)-positivity constraints, we must factor the 4-particle 
expansion coefficients bjkim and djkim into products of 3- 
and 1-particle coefficients bjbkim and bjb* klm which cause 
the cumulant part of the 4- RDM in (*|6|*) to vanish 



b Mmb*b* qst ( 4 A J p k q l ™ 

jklmpqst 



>jklm\ 

pqst 



) = 0. (25) 



The (3,4)-positivity condition, represented by Eq. (|2"Tj) 
and the tensor decomposition of the expansion coeffi- 
cients, is part of a class of (3,4)-conditions that arises 
from all distinct combinations of two 4-particle metric 
matrices that differ from each other in the replacement 
of three second-quantized operators by their adjoints. 

A class of (2,4)-positivity conditions, shown in Table HI 
can be derived from convex combinations of the above 
(3,4)-positivity conditions that cancel the 3-particle op- 
erators, that is the products of six second-quantized op- 
erators. To effect the cancelation, the nonlinearity of the 
expansion coefficients of C must be increased from bjbkim 
to bjCkdie m . Specifically, the C operators in Table U are 
defined as 



= £ b]cld?el 

jklm 



.ajalcii a„ 



(26) 



, l ,, ( , ,,. ,. ., x and (/, and 



where a" and 6" are a] and b* if u 
if u = o. Each of the eight (2,4)-positivity conditions in 
Table U generates an additional condition by switching 
all x's and o's in accordance with particle-hole duality, 
the symmetry between particles and holes. The (2,4)- 
conditions become the diagonal A^-representability con- 
ditions @, UK-Ill] when b, c, d, and e are restricted to 
be unit vectors; they are more general than the unitarily 
invariant diagonal conditions because these four vectors 
are not required to be orthogonal. These (2,4)-positivity 
conditions are only representative of the process by which 
complete conditions can be constructed from the solu- 
tion of the A^-representability problem presented in this 



Letter. Additional (2,4)-conditions in this class can be 
generated from reordering creation and annihilation op- 
erators in the conditions of Table U and other extreme 
(2,4)-conditions can be constructed from lifting the (2,3)- 
conditions. A comprehensive list of (2,4)-positivity con- 
ditions as well as (2,3)-, (2,5)-, and (2,6)-positivity con- 
ditions, which are consistent with the constructive solu- 
tion, will be presented elsewhere [3l|. The (2,5)- and 
(2,6)-conditions include extensions of three and eighteen 
classes of known diagonal conditions, respectively. 

The set Pjy of iV-representability conditions on the 2- 
RDM contains the set C 1 ^- of classical N -representability 
conditions @, lUSl , which ensure that the two-electron 
reduced density function (2-RDF), the diagonal (classi- 
cal) part of the 2-RDM, can be represented by the in- 
tegration of a TV-particle density function. In different 
fields the set C% of A^-representable 2-RDF has been 
given different names: cut polytope [32j in combina- 
torial optimization and the correlation (or Boole) poly- 



tope [33, |36| in the study of 0-1 programming or Bell's 
inequalities. The set C^-, previously characterized, has 
important applications in global optimization including 
the search for the global energy minima of molecular 
clusters [34|], the study of classical fluids [HI], the max- 
cut problem in circuit design and spin glasses (32j, lat- 
tice holes in the geometry of numbers, pair density (2- 
RDF) functional theory [33j, and the investigation of 
generalized Bell's inequalities (36|. The characterization 
of the set of A^-representable 2-RDMs represents a 
significant generalization of the solution of the classical 
A^-representability problem (the Boole 0-1 programming 
problem). In addition to its potentially significant ap- 
plications to the study of correlation in many-fermion 
quantum systems, knowledge of the set Pjy may have 
important applications to "quantum" analogues of prob- 
lems in circuit design and the geometry of numbers. 

The complete set of A^-representability conditions 
firmly solidifies 2-RDM theory as a fundamental theory 
of many-body quantum mechanics with two-particle in- 
teractions. Rigorous lower bounds to the ground-state 
energy of strongly correlated quantum systems can be 
computed and improved in polynomial time from sub- 
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sets of the complete iV-representability conditions [20( 
(Minimizing the energy with a fully 7V-representable 2- 
RDM is a non-deterministic polynomial (NP) complete 
problem because Cf^ C P$r with optimization in 
known to be NP-complete [32j). The present result raises 
challenges and opportunities for future research that in- 
clude (i) implementing the higher iV-representability con- 
ditions which are not in the form of traditional semidefi- 
nite programming 14, 15, 2^], and (ii) determining which 
of the new conditions are most appropriate for differ- 
ent problems in many-particle chemistry and physics. 
Beyond their potential computational applications, the 
complete iV-representability conditions for fermionic den- 
sity matrices provide new fundamental insight into many- 
electron quantum mechanics including the identification 
and measurement of correlation and entanglement. 
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